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The small angular scale fluctuations of the (on large scale) isotropic gamma-ray background (IGRB) carry information 
about the presence of unresolved source classes. A guaranteed contribution to the IGRB is expected from the unresolved 
gamma-ray AGN while other extragalactic sources, Galactic gamma-ray source populations and dark matter Galactic 
and extragalactic structures (and sub-structures) are candidate contributors. 

The IGRB was measured with unprecedented precision by the Large Area Telescope (LAT) on-board of the Fermi 
gamma-ray observatory, and these data were used for measuring the IGRB angular power spectrum (APS). Detailed 
Monte Carlo simulations of Fermi-LAT all-sky observations were performed to provide a reference against which to 
compare the results obtained for the real data set. The Monte Carlo simulations are also a method for performing those 
detailed studies of the APS contributions of single source populations, which are required in order to identify the actual 
IGRB contributors. 

We present preliminary results of an anisotropy search in the IGRB. At angular scales <2° (e.g. above multipole 
155), angular power above the photon noise level is detected, at energies between 1 and 10 GeV in each energy bin, with 
statistical significance between 7.2 and 4.1a. The energy not dependence of the fluctuation anisotropy is pointing to 
the presence of one or more unclustered source populations, while the energy dependence of the intensity anisotropy is 
consistent with source populations having average photon index V = 2.40±0.07. 

Keywords: 



1. IGRB Anisotropy m 

20 

In all-sky high-energy gamma-ray observations an in- 
tense diffuse emission originating from the Milky Way is 22 
visible. Above 30 MeV the large majority of this emission 23 
is produced by cosmic-ray (CR) ions interacting with the 24 
Galactic interstellar gas via neutral pion production and 25 
decay, and inverse Compton (IC) scattering of interstellar 26 
radiation fields photons off CR electrons. A fainter, almost 27 
isotropic on large angular scales, gamma-ray background 28 
(IGRB) has also been detected, first by the SAS-2 mis- 29 
sion [l| . Later the IGRB energy spectrum was measured 30 
with good accuracy by the Energetic Gamma-Ray Experi- 
ment Telescope Q , 0] on-board the Compton Observatory 32 
and recently with the LAT detector on-board the Fermi 33 
gamma-ray observatory Q ■ 34 

A considerable part of the IGRB has likely extragalac- 35 
tic origin (EGB) and carries information about the non- 36 
thermal phenomena in the universe. The IGRB has a guar- 37 

38 
39 
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anted contribution from the known extragalactic gamma- 
ray sources. AGN represent the largest source population 
above 30 MeV detected by EGRET [| and the LAT de- 
tector Therefore, undetected AGN (those below 
the current detection threshold) are the most likely can- 
didates for the origin of the bulk of the EGB emission. 
The estimates of the EGB fraction originating from AGN 
vary from 20 and «100% depending on the energy and the 
model (see 0] and references therein). Another main ex- 
tragalactic candidate is star-burst galaxies [9]. A fraction 
of the IGRB might also originate from the sum of Galactic 
nearby sources such as millisecond pulsars (MSPs) [lbl ]. 

The IGRB anisotropy study provides us also with a 
method for the indirect search for dark matter (DM) with 
gamma rays, which is very complementary to the study 
of the main DM overdensities, such as the Galactic center 
[11| and dwarf spheroidal galaxies [12j and also comple- 
mentary to the search of lines [l3| or other spectral fea- 
tures from large areas of the sky. In fact the sub-structures 
of the Milky Way dark matter halo and the dark matter 
halos in the local universe might produce an imprint in the 
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IGRB. Dark matter overdensities might shine in gamma 95 
rays both in the case of pair-annihilating DM particles 96 
(with the resulting gamma-ray flux proportional to the 97 
square of the dark matter density) and pseudo-stable de- 
caying ones (with flux proportional to the density). The 98 
largest substructures of the Galactic halo might be indi- 99 
vidually detectable while the majority of them are likely to 100 
be under detection threshold but still be able to contribute 101 
to the IGRB [3] and also induce small scale anisotropics 102 
0. " 103 

Information on the IGRB origins is carried by its small 104 
scale angular fluctuations. If a diffuse emission originates 105 
from unresolved source populations then small angular 106 
scale fluctuations will arise because of the variation of 07 
number density of the sources in different sky directions. 108 
These fluctuations are a feature of the source distribution 109 
and will persist also in the limit of infinite statistics, con- 110 
trary to Poisson fluctuations (photon noise) which decre- 111 
ses with increasing event statistics. It is therefore possi- n2 
ble to discriminate between sources induced anisotropics 

L 113 

and photon noise, if one is provided with a sufficient data iu 
statistics. One method for the study of the IGRB anisotropy 
is the APS measurement. In theoretical studies of the 

116 

IGRB anisotropy the following contributors were consid- 
ered: (a) blazars 17|. ll8i.[l9l ; (b) star- forming galaxies 117 
HJ; (c) Galactic MSPs|U(d) annihilatingor decaying 118 
dark matter in Galactic sub-halos | 2 3 . | 2 3| . p3 | : (e) dark 119 
matter in extragalactic structures [16j , [17| , [19( , [24[ , [25j and 120 



28], [2£ 



In these studies it is also shown that 121 

69 intrinsic clustering of many populations candidates has 122 

70 a sub-dominant effect on the angular power spectrum in 123 

71 multipole range between 100-500. 124 

72 Here we report on a search for anisotropies in the IGRB, 125 

73 performed with the data taken with the LAT detector on- 126 

74 board of the Fermi satellite. 127 

128 

75 2. The LAT and the Data Analysis 129 

130 

76 The Fermi Large Area Telescope has a wide field of 

77 view (2.4 sr) and a large effective area («8000 cm 2 for i 

78 normally- incident photons above 1 GeV). LAT is a pair- i 

79 conversion telescope with a modular structure made of 4 x 4 

so towers. Each tower is composed by: (i) a precision silicon™ 

si tracker (18 planes of Si-strip detectors coupled with W135 

82 conversion planes, with a total of 1.5 Xo for the normal in-136 

83 cidence. The Si tracker is made of a first thinner segment, 137 

84 called front, with a better angular resolution and a seconds 

85 thicker one, called back) ; (ii) a Csl homogeneous calorime-139 

86 ter (8.5 Xo for the normal incidence). The pair-conversionwo 

87 telescope is covered with an anti-coincidence detector thatui 
as allows for rejection of charged particle events. Full detailsi42 

89 of the instrument, including technical information about 

90 the detector, on-board and ground data processing, and 143 

91 mission-oriented support, are given in |3lj . 144 

92 Data taken during the first 56.6Ms of observation («22 145 

93 months) were used. The experimental data and Simula- 146 

94 tions were analyzed with the LAT analysis software Sci- 147 



ence Tools version v9rl5p4 with P6_V3 LAT instrument 
response functions (IRFs). 

The main analysis steps were: 

• Data Preparation. By means of the gtselect tool: (i) 
events of diffuse class were selected and with energy 
between 1 and 50 GeV; (ii) data with zenith angle ex- 
ceeding 105° were rejected to reduce Earth gamma- 
ray albedo contamination; (iii) events converted in 
the front and back tracker segments have also been 
divided in order to be analyzed separately. Periods in 
which with LAT was in the South Atlantic Anomaly, 
not in survey mode or with rocking angle exceeding 
52° were discarded with the tool gtmktime. The in- 
tegrated livetime was calculated using the gtltcube 
tool (photon injection step size cos(0) = 0.025, pixel 
size of 0.125° corresponding to a HEALPix [30| map 
with Nside = 512 resolution) 

• Exposure. Exposure maps were calculated using the 
gtexpeube tool with the same pixel size as for gtltcube, 
for 42 logarithmic energy bins spanning 1.04 to 50 GeV, 
in order to have a good knowledge of the energy de- 
pendence of the exposure. 

The event shuffling technique is an alternative method 
for the exposure calculation that does not rely on 
the Monte Carlo based calculation of the exposure 
implemented in the Science Tools. It was used for 
cross-checking possible exposure systematic errors. 
With this method the arrival directions of pairs of 
detected events are swapped in the instrument co- 
ordinates, this produces anexposure map wih arbi- 
trary normalization. The same technique has also 
been used to search for anisotropy in the CR elec- 
trons arrival directions, measured with the LAT, in 
& 

• Intensity Maps. Counts maps were built with the 
selected data. Both the photon counts and exposure 
maps were converted into HEALPix-format maps with 
Nside = 512, and HEALPix gamma-ray intensity 
maps in four energy bins were obtained. 

• Map Masking. Regions of the sky heavily contami- 
nated by Galactic diffuse emission were excluded by 
masking Galactic latitudes \b\ < 30°, and masking 
sources in the Fermi 11-month catalog [ff within a 
2° angular radius. In this study we focused on multi- 
poles Z > 150 (corresponding to angular scales < 2°), 
since lower multipoles (corresponding to correlations 
over larger angular scales) are likely more contami- 
nated by Galactic diffuse emission. 

• Intensity and Fluctuation APS. The angular power 
spectra were calculated for the intensity map I(V0, 
where ip is the sky direction. The angular power 
spectrum is given by the coefficients C/ = < |a/ m | 2 > 
with the ai m determined by expanding the map in 
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148 spherical harmonics. The intensity APS indicates 

149 the dimensionful amplitude of the anisotropy and 

150 can be compared with predictions for source classes 

151 whose collective intensity is known or assumed. A 

152 dimensionless fluctuation APS has also been calcu- 

153 lated by dividing the intensity APS spectrum C; of a 

154 map by the mean sky intensity (outside of the mask) 

155 squared <I> 2 . In the case of the shuffling technique 

156 only the intensity fluctuation APS can be obtained. 

157 The angular power spectra of the maps were calcu- 

158 lated using the HEALPix package [30(. The mea- 

159 sured APS were corrected for the power suppres- 

160 sion due to the beam and pixel window functions, 

161 and an approximate correction, valid at multipoles 

162 I >100, was applied to account for the reduction in 

163 angular power due to masking. For each energy bin, 

164 the APS of the maps of front- and back-converting 

165 events were calculated separately and then combined 

166 by weighted average. 

16? • Validation Studies. Careful checks of the procedures 

168 and the related parameters were performed as well as 

169 tests for assessing the origin of the measured angular 

170 power. The full analysis was applied to a simulated 

171 point source population, in order to compare the 

172 APS determined with the analysis chain and the one 

173 which could be analytically calculated. The depen- 

174 dence of the results on the instrument response func- 

175 tions (IRFs), as also on the masking have also been 

176 studied. For the latter in particular both the latitude 

177 cut with which the Galactic plane was excluded and 

178 the radius of the circle with which each know source 

179 was masked, were varied. In order to further evaluate 
lso the effect of residual Galactic diffuse emission on the 
lsi APS the analysis was repeated on the data after the 

182 subtraction of a model of the Galactic diffuse emis- 

183 sion (Galactic Foreground Cleaning). The details of 

184 these studies will be given in the final publication of 

185 the anisotropy study. 



3. Simulated Models 

Detailed Monte Carlo simulations of Fermi-LAT all- 
sky observations were performed. The purpose was to 
have a reference model to be compared with the real data 
set. The status-of-art of the LAT modelling was used by 
means of the gtobssim tool. This code required as an in- 
put a detailed spatial and spectral model of the emission 
to be simulated. Details of real LAT observations, such 
as spacecraft pointing and live-time history, can be also 
included in the simulations. The gtobssim tool generates 
simulated photon events. The P6_V3_DIFFUSE IRFs and 
the actual spacecraft pointing and live-time history match- 
ing the observational time interval of the data were used 
to generate the simulated data sets. Two models of the 
gamma-ray sky were simulated. Each model is the sum of 
three components: 
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Figure 1: PRELIMINARY. Fluctuation APS of the IGRB minus 
photon noise, in four energy bins. For isotropic emission this differ- 
ence would be consistent with zero. The large angular power at I < 
155 is likely originating from Galactic diffuse emission. For I > 155 
the measured power excess is approximately constant in multipole, 
suggesting that it originates from one or more unclustered source 
populations. The APS obtained with the event shuffling technique 
are also reported. The two different methods provide consistent re- 
sults and show that possible inaccuracies in the exposure calculation 
have a negligible impact on these results. 
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Figure 2: PRELIMINARY. Intensity APS of the IGRB minus photon 
noise, in four energy bins. APS of the experimental data, simulated 
default model and high resolution models are reported. The mea- 
sured power above 1 ?s 155 is not found in either of the two models. 
Points from different data sets are offset slightly in multipole for 
clarity. 
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Figure 3: PRELIMINARY. Intensity APS of the IGRB minus photon 
noise, in four energy bins. The APS of the single model components 
are reported. The large power in the models at I < 155 is due to 
the GAL component. The isotropic component (ISO) provides APS 
compatible with zero as expected for isotropic emission; likewise, 
the source component (CAT) provides no contribution because all 
sources were effectively masked. Total Model points are offset slightly 
in multipole for clarity. 
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1. GAL - a model of the Galactic diffuse emission; 25s 

2. CAT - the 1 451 sources in the first Fermi-LAT sourcezss 
catalog (1FGL) @; 

3. ISO - an isotropic background. 261 

262 

The same CAT and ISO components are included in 
both models. These differ only for the model for the GAL 
component in use. GAL describes both the spatial dis- 263 
tribution and the energy spectrum of the Galactic dif- 
fuse emission. The GAL component for the reference sky 265 
model used in this analysis (hereafter, Model) is the recom- 2 
mended Galactic diffuse model for LAT data analysis, gll 267 
iem_v02.fit 32j, which has an angular resolution of 0.5°. 268 
This model was used to obtain the 1FGL catalog; a de- 

r 1 °' 269 

tailed description can be found in [33j. 270 
An higher-resolution model (Hi-Res Model) was simu- 2n 
lated for comparison, so that was possible to test impact of 272 
smaller details in the Galactic diffuse emission. This model 273 
(ring 2Imonth vf.fit) is internal to the LAT collaboration,^ 
and was built using the same method as gll iem_v02.fit, 
but differs primarily in the following ways: (i) this model 27 
was constructed using 21 months of Fermi -LAT observa- 

277 

tions, while gll iem_v02.fit was based on 9 months of data; 278 
(ii) the grid angular resolution for this model is 0.125° , in 2?g 
order fully exploit the angular resolution of the CO maps 



351 ] used to build it; and (iii) additional large-scale struc- 



tures, such as the Fermi bubbles [36], are included in the 

1 I 'J ' 282 

model through the use of simple templates. 

The single power law spectrum is assumed for all the 2g4 
sources in CAT and the locations, average integral fluxes, 28 
and photon spectral indices are as reported in the lFGL 2g6 
catalog. All 1451 sources were included in the simulation. 2g? 
ISO represents the sum of the Fermi-measured IGRB and 
an additional isotropic component presumably due to un- 2g9 
rejected charged particles; for this component the spec- 290 
trum template isotropic iem_v02.txt was used. For both 29i 
the Model and the Hi- Res Model , the sum of the three sim- 

1 292 

ulated components results in a description of the gamma- 
ray sky that closely approximates the angular-dependent 294 
flux and energy spectrum of the all-sky emission measured 295 
by the Fermi -LAT. Although the simulated models may 
not accurately reproduce some large-scale structures, e.g., 
Loop I [37J and the Fermi bubbles [jjfj, these features are 29g 
not expected to induce anisotropics on the small angular 
scales on which we focus in this work. 300 
The simulated models were processed and their angular 3oi 
power spectra calculated using the same analysis pipeline 302 
as the data. In Figure 2 the APS of the data and mod- 

& 303 

els are shown. The contributions to APS of the individual 

304 

components of the default model are shown in Figure 3. At 
all energies the only component of the models contribut- 306 
ing significantly to / < 155 is the Galactic diffuse emission. 
The contribution from the isotropic component is negligi- 308 
ble, since this component is isotropic by construction and 3og 
thus, after the photon noise is subtracted, it should only 3io 
contribute to the monopole (^=0) term. The source catalog^ 
component contributes zero power at all energies and mul- 



tipoles because all these sources were effectively masked, 
giving rise to a negligible residual effect. We remark that 
in general the APS of distinct components are not linearly 
additive due to contributions from cross-correlations be- 
tween the components. 

4. Results 

Significant (>3<r CL in each energy bin) angular power 
above the photon noise level (see Figure 1) is detected in 
the data at multipoles 155 < I < 504 for energies below 
10 GeV, and at lower significance in the 10-50 GeV energy 
bin. The sensitivity of the measurement at high energies 
is limited by the small statistics. The measured angular 
power is consistent with a constant value within each en- 
ergy bin, which suggests that it originates from one or 
more unclustcrcd populations of point sources. The fluc- 
tuation angular power detected in this analysis falls below 
the level predicted for some models of blazars, MSPs, and 
Galactic and extragalactic DM structures, and so the mea- 
sured amplitude of the fluctuation angular power can limit 
the contribution to the total IGRB intensity of each source 
class. The measured fluctuation angular power is almost 
independent of energy, and so it might originate from a 
single dominant source class, although a mild energy de- 
pendence cannot be excluded. 

The energy dependence of the intensity angular power 
of the data is well-described by that arising from a single 
source class with a power-law energy spectrum with pho- 
ton index T= 2.40± 0.07. This value is compatible with 
the mean intrinsic spectral index for blazars as determined 
from recent Fermi-LAT measurements. 

The study of the currently available Fermi-LAT data 
has provided us with the first detection of small scale- 
angular power in the IGRB. Future analyses on larger 
data sets could allow anisotropy searches to be extended 
to higher energies, and will likely be more sensitive to fea- 
tures in the anisotropy energy spectrum. 

Detailed studies of population models are required in 
order to: (i) identify the source classes which are actu- 
ally contributing to the IGRB anisotropy and (ii) provide 
upper-limits to the contribution of other candidate source 
populations. A possible approach is the study of the APS 
induced by source population models and the following 
Monte Carlo simulation of the model emission, with the 
best possible characterization of the instrument response. 

For an example, in the DM case a model of all its dis- 
tribution structures and sub-structure (Galactic and in the 
local Universe) which can provide contribution should be 
considered [38j |. and furthermore a given particle physics 
model should be assumed, or at least the DM particle mass 
and the annihilation cross-section for self-annihilating DM, 
or the mean life-time for pseudo-stable decaying DM. Then 
the APS induced by these DM models can be obtained 
with the same procedure which was used for the simulated 
Model and Hi- Res Model of section 3. 
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